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Abstract 

We present results of a systematic study of the c axis lattice dynamics in single layer Bi2Sr2Cu06 
(Bi2201), bilayer Bi2Sr2CaCu208 (Bi2212) and trilayer Bi2Sr2Ca2Cu30io (Bi2223) cuprate super- 
conductors. Our study is based on both experimental data obtained by spectral cUipsomctry on 
single crystals and theoretical calculations. The calculations are carried out within the framework 
of a classical shell model, which includes long-range Coulomb interactions and short-range interac- 
tions of the Buckingham form in a system of polarizable ions. Using the same set of the shell model 
parameters for Bi2201, Bi2212 and Bi2223, we calculate the frequencies of the Brillouin-zone center 
phonon modes of Pl2u symmetry and STiggcst the phonon mode eigenvector patterns. We achieve 
good agreement between the calculated A2u eigenfrequencies and the experimental values of the c 
axis TO phonon frequencies which allows us to make a reliable phonon mode assignment for all 
three Bi-based cuprate superconductors. We also present the results of our shell model calculations 
for the F-point Ai^ symmetry modes in Bi2201, Bi2212 and Bi2223 and suggest an assignment that 
is based on the published experimental Raman spectra. The superconductivity-induced phonon 
anomalies recently observed in the c axis infrared and resonant Raman scattering spectra in trilayer 
Bi2223 are consistently explained with the suggested assignment. 

PACS numbers: 74.72.Hs, 74.25Kc, 78.20.Bh, 63.20.Dj 



*Also at the Institute for Solid State Physics, Russian Academy of Sciences, Chernogolovka, Moscow distr., 
142432 Russia 
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I. INTRODUCTION 



The unusual c axis charge transport in the cuprate high-Tc superconductors is still a 
subject of a controversial discussion. Of particular interest is a sizable absorption peak that 
develops below Tc in the far-infrared (FIR) range, accompanied with the strongly anomalous 
temperature dependence of some of the c-polarized IR phonon modes, in compounds that 
contain more than one CUO2 plane per unit cell. Recently a reasonable description of these 
effects has been obtained with the so-called Josephson-superlattice model 
Here it is assumed that the individual CUO2 planes are only weakly coupled by Josephson 
currents in the superconducting state. For bilayer compounds like YBa2Cu307_5 (Y123) 
and Bi2Sr2CaCu208 (Bi2212) this results in two kinds of Josephson junctions with different 
longitudinal plasma frequencies, and their out-of-phase oscillation gives rise to a transverse 
Josephson plasma resonance, which has been assigned to the absorption peak that develops 
below Tf,; the Josephson currents lead to a modification of the dynamical local electric 
field. a.d *us a. at Kea.. of.. ob.e„ed phono. a„o.aUe. j QBQ. The t.laye. 
compounds provide an interesting testing ground for the various models that have been 
put forward to explain the unusual c axis charge dynamics and the anomalous temperature 
dependence of some of the IR-active c-polarized phonon modes such as observed in the bilayer 
compounds. Indeed, very recent experimental study of the c axis optical conductivity of a 
trilayer Bi2Sr2Ca2Cu30io (Bi2223) provides clear evidence that the transverse Josephson 
plasma resonance is a universal feature of the multilayer high- TV cuprate compounds Ql- 
It has been shown that the Josephson-superlattice model [2, 0,0,01 allows qualitative 
description of the transverse Josephson plasma resonance and related phonon anomalies in 
Bi2223, which suggest that the Josephson currents lead to a strong variation of the dynamical 
local electric field even between the inner and the outer CUO2 planes of a trilayer. 

The assignment of the IR-active c-polarized phonon modes and their polarization dia- 
grams are of great importance for understanding these phenomena. While the properties 
of the IR-active phonon modes are indeed well known for the Y123 system 



the properties 

Qy,y,iiii, 



the situation is not that clear for the bilayer high- Tp comp ounds of the ty pe A^S r^CaCu^Os 



(A2212; A = Tl, Bi). Despite of experimental p, |6|, 1121, Il3| and theoretical [ij, ll5|, |16| efforts 



/• pe A^a r^ua 



the assignment of the c axis IR phonon modes in Bi2212 and T12212 remains contradictory, 
as discussed in The situation is even less satisfactory for the trilayer high-Tc compounds 
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A2Sr2Ca2Cu30io (A2223) where the lack of single crystals of sufficient size and quality has 
inhibited detailed experimental investigations. Shell model calculations have been carried 
out for T12223 17f and IR measurements have been performed on ceramic samples of T12223 

n n 

UM and more recently on oriented ceramics of Bi2223 [19|. The situation has changed only 
very recently when some of the present authors managed to grow high quality single crystals 
of the trilayer Bi-compound Bi2Sr2Ca2Cu30io that are suitable for IR measurements with 
the light polarization along the c axis direction 2o| . 

This manuscript reports a systematic investigation of the c axis lattice dynamics of the 
bismuth-based compounds Bi2Sr2Ca„_iCu„02n+4+d that contain either n = 1, 2, or 3 Cu02 
layers per unit cell. Our approach combines theoretical shell model based calculations with 
ellipsometric measurements of the c axis dielectric response on well characterized single 
crystals of high quality [2^ 0, l^]. From the ellipsometry experiments we obtain the 
frequencies of the c axis polarized IR-active phonon modes as well as the dielectric constants. 
On the theoretical side we model the Bi-based compounds with n = 1, 2, and 3 in the 
14/mmm space group using the shell model approach suggesting that screening due to free 
charge carriers should not have major effects [2^. This approach has been successfully 



a| 

L 



ied to study lattice dynamics and various physical properties of high-T^ superconductors 
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2J,|25|. The calculations have been performed using LP code Q- 

11 

model parameters have been derived by fitting to the crystal structures in equilibrium as well 
as to experimental values of the dielectric constants and frequencies of the transverse optical 
(TO) phonon modes. Using the same set of shell model parameter for the n = 1, 2, and 
3 compounds we calculate the eigenfrequencies of the Brillouin-zone center phonon modes 
of A2M symmetry and derive the phonon mode eigenvector patterns. Our results for the 
bilayer compound Bi2212 are consistent with those of the earlier shell model calculations 
in the I^/mmm space group by Prade et al. Q|. We make an assignment of the c axis 
polarized phonon modes observed in the ellipsometric spectra of single, bi-, and trilayer 
Bi-based compounds by comparing with the calculated eigenfrequencies of A2M symmetry 
modes. In addition, we also show the displacement patterns and calculated eigenfrequencies 
for the Raman-active phonon modes of Aig symmetry according to the results of the present 
shell model calculations. 
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II. CRYSTAL STRUCTURE AND T-POINT PHONONS IN BiaSrsCan iCun02n+4+5 
WITH N = 1, 2 AND 3 

As a first approximation, a body-centered-tetragonal structure {M/mmm) has been most 
frequently used in the interpretation of the IR and Raman spectra of Tl- and Bi-based 
compounds. According to the x-ray data 28], the three Bi-based cuprates of general 
formula Bi2Sr2Ca„_iCu„02n+4+(5 with n = 1, 2 and 3 have similar structures of the tetrag- 
onal I4/mmm space group [atetr = 3.814 A), with two formula units (two primitive unit 
cells) in the crystallographic cell, as illustrated in Figure 1. Primitive unit cells of single, 
bi- and trilayer Bi-based cuprates differ only in the number (n-1) of Cu02-Ca-Cu02 slabs 
packed along the c axis; upon insertion of one and two slabs the c axis parameter of the 



crystallographic cells increases from 24.6 to 30.6, and 37.1 A 28|. The Wyckoff positions 
of the atoms and their site symmetries in tetragonal lA/mmm (D^^) primitive unit cells of 
Bi2Sr2Cu06 (Bi2201), Bi2Sr2CaCu208 (Bi2212), and Bi2Sr2Ca2Cu30io (Bi2223) are listed 
in Table 1, with the notations of the atoms given in Figure 2. 

The irreducible representations corresponding to various atomic sites in the three Bi- 
based compounds that follow from the character tables of the point groups are presented 
in the right column of Table 1. For the IR- active modes at the F-point of the Brillouin 
zone the displacement symmetries are A2„ for a displacement in the z direction and 
for displacements in the x or y directions, whereas those for the Raman-active F modes 
are Ai^, Bi^, and E^,, respectively. Due to full point-group symmetry of the Cul (Bi2201), 
Cu2 (Bi2223) and Ca (Bi2212) atoms with the (la) site position, no Raman-active modes of 
these atoms are allowed, and the three displacement vectors yield three modes of and E^^ 
(doubly degenerate) character. The oxygen atoms in the mirror cuprate planes, 01 in single 
layer and 04 in trilayer Bi-compounds, of site symmetry, yield eigenmodes + B2m + 
2Eu. The in-plane oxygen 01 with C2i, site symmetry in the bilayer and trilayer compounds 
generates eigenmodes A2M + Ai^ + Big + + 2Ep + 2E„, among them the mode 
is silent and E^, is doubly degenerated. The total numbers of modes grouped according to 
their optical activity in the three Bi-based compounds are summarized in Table 2. 

It has been shown that the actual structure of Bi-based compounds is better represented 



by orthorhombic Cccm (D^^) |29|, l3Cll I31L l32|| or Ccc2 (Ca^)) [331 134| space groups, most of 



the structural studies verify that the orthorhombic Cccm structure is correct. To follow the 
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TABLE I: The irreducible representations for the atoms in tetragonal I^/mmm (04^) Bi2Sr2Cu06 
(Bi2201), Bi2Sr2CaCu208 (Bi2212) and Bi2Sr2Ca2Cu30io (Bi2223). 



Comp. 


Atom 






Wvc 
not. 


Site 
svm 


Irrednri V)le 


rPDresentation 




111 








Ft., 


A2u 


J- F 




Bi, Sr, 


02, 03 




f2e) 


04f 








01 








Do;, 






Bi2212 


Ca 






fla) 


'-'Ah 








Bi, br, 


Oul, 02, 


03 


(2e) 




5(A2„ + Alg 


+ Eg + E„) 




01 






(4g) 


2v 


A2U + Alg + Blc, 


+ B2u + 2Eg + 2E„ 


Bi2223 


Cu2 






(la) 


^Ah 


^2u 


+ E« 




04 






(2c) 


^2h 




+ B2u + 2Eu 




Bi, Sr, 


Ca, Cul, 


02, 03 


(2e) 




6(A2„ + Alg 


+ Eg + E„) 




01 






(4g) 


2v 


A2U + Alg + Big 


+ B2u + 2Eg + 2Eu 



TABLE II: Modes classification from the irreducible representations for Bi2201, Bi2212 and Bi2223 
in tetragonal I4/mmm (D\1) space group. 

Comp. Tin Tn 

aman ^acoustic ^silent 

Bi2201 5A2„ + 6E„ 4Aig + 4E3 A2„ + E„ B2n 

Bi2212 6A2„ + 7E„ 6Aig + 7Eg + Big A2„ + E„ Bg^ 

Bi2223 8A2„ + 10E„ 7Aig + 8Eg + Big A2„ + E„ 2B2„ 



consequences of the lowering symmetry it is convenient to use the Amaa space group, which is 
a nonconventional setting of Cccm. The orthorhombic Amaa structure has two formula units 
per primitive face centered cell, with in-plane lattice parameters a and b of approximately 
atetrV^, and four formula units for crystallographic cell, with the same c axis parameter. By 
way of example, in the right column of Table 3 we present the irreducible representations 
corresponding to various atomic sites in the orthorhombic Amaa Bi4Sr4Cu20i2 (Bi2201) 
and Bi4Sr4Ca2Cu40i6 (Bi2212) according to the structural data of Imai et al. [29], where 
we also show the modes which are silent in the tetragonal symmetry. The total numbers 
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TABLE III: The irreducible representations for the atoms in orthorhombic Amaa (Dife) 
Bi4Sr4Cu20i2 (Bi2201)and Bi4Sr4Ca2Cu40i6 (Bi2212). 

Comp. Atom Wye. Irreducible representations 

not. 

Bi2201 Cul (2e) A„ + 2Bi„ + 2B2„ + Bg^ 

Bi„ + B2u tetragonal silent 
Bi, St, 02, 03 (41) 4(A„ + 2Bi„ + 2B2„ + Bg^ + 2Ag + Bi^ + Bsp + 2B35) 

4(Bi„ + B2„ + Ag + Big + Bag) tetragonal silent 
01 (4h) A„ + Bi„ + 2B2„ + 2B3„ + A^ + Bi^ + 2B2g + 2B3g 

Ag + Big + 2B2g + 2B3g tctragoual silent 
Bi2212 Ca (2e) A„ + 2Bi„ + 2B2„ + B3„ 

Bi„ + B2,/ (.('(ragoual silcut. 
Bi, Sr, Cul, 02, 03 (41) 5(A„ + 2Bi„ + 2B2„ + B3„ + 2Ag + Bi^ + Bsg + 2B3g) 

5(Bi„ + B2« + Ag + Big + Bo^g) tetragoual silent 
011,012 (4h) 2(A„ + Bi„ + 2B2„ + 2B3„ + A^ + Bi^ + 2B2g + 2B3g) 

Biu + 2B2u + 2B3„ + Ag + Big + 2B2p + 2B3g tetragonal silent 



of modes grouped according to their optical activity in the Amaa Bi2201 and Bi2212 are 
summarized in Table 4. In Table 5 we show components of the IR and Raman-active 
normal modes of the atoms with different Wyckoff positions at the T-point of the Brillouin 
zone in the orthorhombic Amaa space group. Prom Tables 1 and 3 one can see that if 
the 14/mmm symmetry breaks down owing to the orthorhombic distortion of the structure, 
the degeneracy is removed. Thus, for each atom with (2e) site position in the 14/mmm 
symmetry (Bi, Sr, 02, 03 and Cul in Bi2212) normal modes A2U, Aig, one E^, and one 
Eg will become nondegenerate and, as a result, split into two; in addition, one new mode 
of Big character will appear in the orthorhombic Amaa structure. For the atoms Cul in 
Bi2201 and Ca in Bi2212 two new IR-active modes of Bi„ and B2U character arc expected 
in the Amaa structure. Due to breaking of the tetragonal symmetry the 01 atom in Bi2201 
will become " Raman- active" , and, as a result, a number of Raman-active modes of A^, Bi^, 
2B2g, and ^B^g character will appear. 
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TABLE IV: Modes classification from the irreducible representations for Bi2201 and Bi2212 in 
orthorhombic Amaa space group. 



Comp. 




^ Raman 


^acoustic 




^silent 


Bi2201 


10Bi„ + 11B2„ +6B3„ 


9Ag + bBig + 6B25 + lOBgg 


Blu + B2m 


+ Bsu 


6A„ 


Bi2212 


13Bi., + 15B2. +9B3. 


l2Ag + 7Blg + m2g + I4B39 


Blu + B2u 


+ Bsu 


8A„ 


TABLE V: Components of normal 


modes of atoms with different Wyckoff positions 


in orthorhombic 


Amaa space group. 










Wye. not. 


Kg Au 


Bsg Big 


Blu 


B23 


B2m 


(2e) 


X 


X 


zy 




yz 


(41) 


zy X 


yz X X 


zy 


X 


yz 


(4h) 


z z 


xy xy z 


z 


yx 


yx 



However, it has been known that the actual structure of Bi-based cuprates is only pseudo- 
orthorhombic due to stacking faults, extra-stoichiometric oxy gen and incommensurate super- 
structural modulations along the b axis in BiO and SrO layers , which result in substantial 
oxygen and cation disorder. These factors can yield additional phonon modes, in particular, 
"disorder-induced" modes. The crystal structure of Bi-based compounds therefore should 
be described in terms of a larger elementary cell and lower crystal symmetry. The origin 
of the superstructural incommensurate modulations of the atoms in BiO and SrO layers is 
not fully understood. The most common assumption is that the modulations are associated 
with extra-stoichiometric oxygen in the bilayer Bi202+<5 that affects the orthorhombicity in 
the Bi-based cuprates. In this situation, whatever approximation either tetragonal or or- 
thorhombic is adopted to make an assignment in the experimental phonon spectra, the first 
problem one has to solve is to distinguish the extra phonon modes - beyond the considered 
symmetry - due to effects of superstructural incommensurate modulations and presence of 
extra-stoichiometric oxygen in the Bi-based compounds. 

Since the deviation from tetragonality in Bi-based compounds is not large we suggest 
to start with the tetragonal M/mmm approximation assuming that modes silent in tetrag- 
onal symmetry will be weak, and the simple tetragonal lA/mmm approximation should 
provide a useful first approximation for the assignment of the most intense IR and Raman- 
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active phonon modes. To explain more details in the experimental spectra due to possible 
orthorhombicity effects we suggest to invoke the group theory analysis based on the or- 
thorhombic Amaa (Dg^) space group (see Table 3). 



III. SAMPLES AND EXPERIMENTAL TECHNIQUE 

The Bi2201, Bi2212 and Bi2223 single crystals have been grown by the travelling solvent 



ffoating zone (TSFZ) method [2^|2l|, |2^. The feed rods with the nominal composition of 
Bi2.iSri.9Ca„_iCu„02n+4+<5 (n = 1, 2 and 3) were prepared by the conventional solid state 
method. The slow growth rates of 0.5 and 0.2 mm/h were used for the Bi2201 and Bi2212 
crystals, respectively. Since the Bi2223 crystallization field is the narrowest among the three 
members of Bi-based family, an extremely slow growth rate of 0.04 mm/h was used. The 
Bi2201, Bi2212 and Bi2223 crystals were grown under applied oxygen pressure p(02), in 
air, and in mixed gas flow of 80% Ar and 20% O2, respectively. Because of the significant 
difference between the lattice parameters of the a and h axis as compared to the one along 
the c axis, the growth rate is highly anisotropic and very slow along c. One typically obtains 
platelike crystals with large ah planes of about 6x2 mm^ but with much smaller dimensions 
along the c axis. The typical thickness is 1.5 2 mm for Bi2201, about 1 mm for Bi2212, 
less than 1 mm for Bi2223. 

The crystals were characterized by energy dispersive X-ray (EDX) analysis. X-ray diffrac- 
tion (XRD), transmission electron microscopy (TEM), magnetic susceptibility and resistivity 
measurements. The as-grown Bi2212 crystals are phase-pure with a pseudotetragonal struc- 
ture, they are almost optimally doped with Tc = 91 K and ATc = 2 K. The as-grown 
Bi2223 crystal contains 95-98 % of the Bi2223 pseudotetragonal phase with a minor fraction 
of layer- intercallated Bi2212. In the as-grown state the Bi2223 crystal was underdoped with 
Tc = 97 K (midpoint) and ATc (10-90% shielding fraction) = 7 K. After a post-annealing 
treatment for ten days in flowing oxygen at 500 °C with a subsequent rapid quenching to 
room temperature, the crystal became nearly optimally doped, with Tc = 107 K and ATc = 3 
K. The undoped Bi2201 crystals grown by the TSFZ method are often nonsuperconducting. 
The superconductivity of the as-grown Bi2201 crystals is critically controlled by the oxygen 
content, therefore the growth was carried out under oxygen pressures p{02) ranging from 1 
to 10 bar. The Tc of the as-grown Bi2201 crystals, determined from magnetic susceptibility 
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measurements varied between 3.5 and 8.5 K upon the Bi/Sr ratio and the oxygen pressure 
p{02) apphed in the course of growing. The structure of the as-grown Bi2201 crystals was 
determined to be pseudotetragonal, the lattice parameters were found to dependent upon 
the Bi/Sr ratio in the crystal composition. A pair of as-grown Bi2212 single crystals (taken 
from the same part of an ingot) was used for oxygen isotope replacement. The selected 
samples were annealed in isotopically enriched oxygen under identical conditions following 
a similar procedure described in Ref. js^. 

For the ellipsometric measurements the crystal surfaces were polished to optical grade. 
The technique of ellipsometry provides significant advantages over conventional reflection 
methods in that (i) it is self-normalizing and does not require reference measurements and 
(ii) €i{u) and e2{uj) are obtained directly without a Kramers-Kronig transformation (see 
Refs. flEl for a description of the technique). The ellipsometric measurements have 
been performed at the U4IR beamline of the National Synchrotron Light Source (NSLS) 
at Brookhaven National Laboratory, USA and at the infrared beamline of the synchrotron 
radiation source ANKA at Karlsruhe Research Center, Germany. Home-built ellipsometers 
attached either to a Nicolet Magna 860 or a "Bruker" IFS 66v/S FT-IR spectrometer were 
used. The high brilliance of the synchrotron light source due to the small beam divergence 
enabled us to perform very accurate ellipsometric measurements in the far-IR range even on 
samples with comparably small ac faces. 



IV. C AXIS IR SPECTRA OF Bi2201, Bi2212, AND Bi2223 SINGLE CRYSTALS 

Figure 3 shows the real parts of the c axis optical conductivity o"i(c(j) and the dielectric 
function ei(ti;) at 150 K for the Bi2212 and Bi2223 single crystals with = 91 K and = 
107 K. It is evident that the optical spectra look rather similar for both compounds. This 
is a consequence of the generic structural features of these two compounds. The electronic 
background is extremely weak and the spectra are dominated by the contributions of several 
IR-active phonon modes. Using the classical dispersion analysis we fit a set of Lorentzian 
oscillators simultaneously to £i{uj) and e2{uj) in the investigated spectral range. To describe 
the weak and rather flat electronic background in a Kramers-Kronig consistent way we have 
also included a sum of a Drude term and several very broad Lorentzian oscillators. As for 
the present purpose this description is rather formal we only limited the half-widths of the 
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TABLE VI: Frequencies uj (in cm~^), oscillator strengths Sj, and damping parameters 7j (in cm^^) 
of the TO phonon bands observed in the c-polarized £2(0;) spectra (T = 150 K) of the single crystals 
Bi2201, Bi2212, and Bi2223. The assignment with the calculated frequencies of {I4:/mmm) 
eigenmodes is given. 



Bi2212 


7j 


A2« 


Bi2223 


7j 


Mu 


Bi2201 


7j 




97 


- 


93 


97 


- 


88 


109 (0.56) 


5.00 


77 


- 


- 


- 


127 (0.10) 


9.1 


145 


- 


- 


- 


168 (0.29) 


8.7 


195 


170 (0.63) 


11.2 


211 


165 (0.45) 


12.3 


175 


210 (0.34) 


12.7 


304 


212 (0.77) 


21.4 


252 


208 (0.65) 


14.5 


- 


282 (0.16) 


12.3 


- 


277 (0.50) 


17.4 


- 


- 


- 


- 


304 (1.13) 


21.7 


343 


305 (1.43) 


27.6 


335 


301 (1.31) 


24.7 


293 


358 (0.33) 


30.7 


403 


360 (0.46) 


30.7 


378 


385 (0.85) 


25.2 


353 








402 (0.16) 


24.0 


516 








462 (0.02) 


10.5 




466 (0.02) 


10.5 










521 (0.22) 


85.3 




512 (0.26) 


81.0 




501 (0.10) 


48.0 




583 (0.65) 


38.2 


643 


582 (0.61) 


41.9 


666 


586 (0.42) 


31.2 


596 


640 (0.03) 


17.8 




640 (0.03) 


17.8 




639 (0.01) 


12.5 




eoo = 4.91 






£00 = 5.77 






eoo = 4.69 







broad Lorenzians to values between 150 and 500 cm~^. In Table 6 we list the frequencies 
Uj, oscillator strengths Sj, and damping parameters 7^ of the TO phonon bands observed in 
the 62(0;) spectra of Bi2212 and Bi2223 single crystals, together with determined values of 
the high-frequency dielectric constants, e^. 

For the Bi2223 we assign the main phonon bands at 97, 127, 170, 212, 305, 360, 402, and 
582 cm^^ to the eight A2„ character modes. These modes are clearly identified especially due 
to the resonance features in ei(cj) as shown in Figure 3(b). Except for the two new intrinsic 
phonon modes in Bi2223 which appear at 127 and 402 cm~^, the modes are observed at 
similar frequencies in both compounds. We believe that the remaining phonon features in 
Figure 3 around 277, 466, and 640 cm~^ are beyond the tetragonal M/mmm approximation. 
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as discussed in Section 2. The feature at 640 cm~^ appears to be a weak satellite of the pro- 
nounced high-frequency mode at 583 cm^^. It is likely related to the effects of characteristic 
incommensurate modulations which induce strong distortions in the SrO and BiO planes. 
Indeed, this feature almost disappears with Pb doping which significantly reduces these su- 



perstructural modulations 0,|38|. Concerning the mode at 277 cm^^, on the low-frequency 
side of the pronounced band at 305 cm~^, we find that its oscillator strength correlates with 
the amount of extra-stoichiometric oxygen that is incorporated within the Bi202+5 bilayer 
during the oxygen annealing treatment By contrast to the 640 cm~^ satellite this feature 
remains persistent with Pb doping j^Sl- As a very weak feature at around 466 cm~^ is quite 
distant from the main phonon features we associate with the character modes, we would 
suggest that it is due to a contribution from some in-plane oxygen-related (02 or 03) 
mode which has a z{c) component (see Tables 3 and 5), or it can be a " disorder- induced" 
defect mode. 

The 150 K spectra of the real parts of the c axis optical conductivity cri{uj) and dielectric 
function ei{uj) of the as-grown nonsuper conducting Bi2201 single crystal, with the stoichio- 
metric composition Bi2Sr2Cu06, are shown in Figure 4. The values of the fitted frequencies 
Uj, oscillator strengths 5*^-, and damping parameters 7j of the TO phonon bands in the 
€2(1-1^) are also listed in Table 6. To the best of our knowledge this is the first report of the 
c-polarized phonon spectra of the Bi2201 single crystal. Also we would like to note that 
the c axis spectra of the Bi2201 reported here exhibit remarkable difference with respect 
to those of the isostructural single layer Tl compound T12201 Six strong TO phonon 
features at around 109, 165, 208, 301, 385, and 586 cm~^ are clearly identified in the ai{uj) 
spectrum of Bi2201, as shown in Figure 4(a), whereas only five character modes are 
expected from the consideration in the ideal I4/mmm space group (see Tables 1 and 2). In 
addition to a weak satellite feature at around 640 cm~^, which we have already discussed in 
the previous paragraph, one "extra" mode appears at low frequencies in the Bi2201. This 
presumably refiects that the deviation from the tetragonality is most pronounced for the 
single layer Bi2201. We expect significant amplitudes in the c-polarized IR spectrum for the 
tetragonal forbidden Bi„ mode of the Bi atoms, which are known to deviate most strongly 
from the tetragonal symmetry. We can also suggest that one new Bi^j mode associated with 
the (2e) site position of the Cul atoms in the Amaa Bi2201, which is tetragonal silent, can 
appear in the orthorhombic symmetry (see Table 3). Also we have found that the spectra of 
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superconducting Bi2201 crystals with Tc varying between 3.5 and 8.5 K depending on the 
Bi/Sr ratio and the oxygen stoichiometry exhibit almost the same phonon features. 

We have furthermore investigated the oxygen isotope effect for ^^O and ^^O substitution 
in Bi2212 as shown in Figure 5. For a mode with predominantly oxygen character the 
expected red shift of the eigenfrequency should be proportional to [m(^^0)/m(^^0)]^/^ - 1, 
i. e. about 5.7%. A lower value of this so-called isotopic frequency shift signals that the 
corresponding phonon mode contains a significant contribution from the heavier (metal) 
ions. As shown in Figure 5(a), a noticeable red shift upon ^^O substitution is observed for 
the four phonons at 210, 304, 358, and 583 cm~^; the estimated values of the shifts are 
4.1, 4.4, 5.3, and 5.5%, respectively. This result confirms the expected trend that the three 
high frequency modes at 304, 358, and 583 cm~^ have predominantly oxygen character. It is 
an interesting finding that the low-frequency mode at 210 cm~^ also exhibits a sizeable red 
shift, which suggests that it contains a significant contribution from the oxygen vibrations. 

Some of the c-polarized IR phonon modes in the bilayer Bi2212 and trilayer Bi2223 
compounds exhibit strongly anomalous temperature dependence at Tc, accompanied with 
the development of a sizable absorption peak below Tc in the FIR range 0,0,0,0] • Figures 
6(a,b) show the real part cri(co') of the c axis optical conductivity of the Bi2212 and Bi2223 
compounds at two different temperatures, above Tc and well below T^. In the normal state 
the spectra exhibit hardly any noticeable changes, except for a sharpening of the phonons 
with decreasing temperature. Right below Tc, however, the spectra change appreciably. This 
is also illustrated in Figure 6(c) which displays the difference ai{T « Tc, cu) — cri(T > Tc, u). 
The most prominent feature is the broad absorption band around 500-550 cm~^ which 
appears below Tc and grows rapidly with decreasing temperature. This absorption band 
has been identified in both Bi2212 and Bi2223 and has been attributed to the transverse 
Josephson-plasma resonance (t-JPR) that is a universal feature of the multilayer high-Tc 
cuprate compounds P, 0, 0]. It is evident from Figure 6 that the formation of the t-JPR 
is associated with an anomalous temperature dependence of the phonon modes at 360, 402 
(specific for trilayer Bi2223) and 582 cm~^ denoted by A, B and C in Figure 6. The principal 
difference between the phonon anomalies observed in the bilayer and trilayer Bi-systems at 
Tc is that in the latter system the changes of the spectral weight of the apical phonon 
mode are much more pronounced and the additional phonon mode exhibiting the opposite 
temperature dependence is observed at 402 cm~^. 
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V. THEORETICAL CALCULATIONS 



A. Shell model approximation and potential parameters 

To estimate the frequencies at which the various eigenmodes appear and, consequently, 
to make the assignment in the experimental phonon spectra (at T = 150 K), we have per- 
formed lattice-dynamical calculations using the shell model approximation. We performed 
the calculations for the three Bi-based structures Bi2201, Bi2212, and Bi2223 within the 
ideal tetragonal M/mmm (D^^) space gro up , using the lattice constants and internal z 



structural parameters reported in Refs. j27 



28l |. Screening has not been taken into account 



in this modeling. But screening should not noticeably affect the IR phonons along the poorly 
conducting c axis direction which we investigate in this work. Moreover, according to the 
study of the phonon screening of the in-plane IR-active modes in high-T^ superconductors 
the majority of the lattice modes have oscillator strength similar to those found in the un- 
doped insulating materials 39j. This surprising result is explained by the inhomogeneous 
charge distribution in the Cu02 planes due to which the ionic cores are no longer effectively 
screened, and the carriers will not be able to effectively screen the induced dipole moment 
Sgj ]. Therefore, we apply the shell model approach suggesting that that screening due to 
free charge carriers should not have major effects [2^. Earlier the shell model approach 
has been successfully applied to study lattice dynamics and various physical properties of 



high-Tr superconductors 
code 



25| . Our calculations are performed using the GULP 
26| . In the context of the shell model, the lattice is considered as an assembly of 
polarizable ions, represented by massive point cores and massless shells coupled by isotropic 
harmonic forces. The interaction includes contributions from Coulomb and short-range in- 
teractions. Short-range interactions between relatively small cations are usually ignored. 
The core and shell charges and the spring constant of each ion are parameters of the model. 
The short-range potentials used for the shell-shell interactions are of the Buckingham form 

Vij = Aijexp{-r/pij) - Cij/r^. (1) 

The Buckingham potential parameters as well as the atoms' core and shell charges and the 
relevant spring constants were fitted in this work. To obtain a self-consistent set of the shell 
model parameters acceptable for a family of Bi-based high-T^ cuprates, we performed fits 
simultaneously to the structures of Bi2201, Bi2212, and Bi2223. As the starting values for 
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the shell model parameters - pair potentials, net ionic charges, shell and core charges and 



relevant spring constants - we used the ones obtained by Baetzold j2J] for a bilayer cuprate 
YBa2Cu307„5 system. We have two reasons to keep the net ionic charges of the atoms in 
the cuprate planes of +2| e | for Cul and Cu2, and -2| e | for 01 and 04 . First, it is 
in line with the approach of Baetzold for Y123 2J| and second, this is required if we want 
to describe single, bi- and trilayer Bi-systems by the same set of the potential parameters. 
Otherwise, by dividing the non-integer charges between the cuprate planes we would have 
different pair potential parameters in each case. Following the approach of Baetzold, we 
also keep the same effective charges of -1.67| e | for 02 and 03. We take the formal ionic 
charge +2\ e \ for Sr, and to satisfy the electroneutrality conditions we have to assume an 
effective charge +2.33| e | for Bi. The formal charges here describe the dipole moment per 
unit displacement, rather than the charge within some region of space. It was noted by 
Catlow and Stoneham that even the integer charge-distribution method can be successfully 
employed to accurately calculate physical properties of solids where ionicity may vary over 
a range of values [l^. We assume the shell charge of -3.2576 1 e | for all oxygen ions, as in 
[2^ . The Sr, Bi, Ca and O ions are treated as polarizable, while the Cu ions are treated as 
unpolarizable. The fitting procedure keeps 0-0 (shell-shell) interaction as those of typical 



oxides 
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42| . The structures were allowed to relax to equilibrium conditions under the 
symmetry {M/mmm) restrictions, the cores and shells were allowed to move separately. The 
forces on the atoms were minimized by varying the potential parameters, the eigenvalues 
for the first three modes for all three compounds had approximately zero values, indicating 
a stable structure in each case. Normally a good fit requires some information on physical 
properties - elastic, dielectric, and piezoelectric (where applicable) constants - as well as the 
structure. In general, vibrational frequencies contain far more information than any of the 
above. However, in fitting the frequencies a reasonable implication for the assignment has 
to be suggested, at least for some key modes. In the present work a final fit was done to the 
available experimental data on the bilayer Bi2212 compound. The criterion of a successful 

n 

fit was good agreement with the structural parameters [27], provided that approximately 
zero forces are applied for all atoms, with the c axis static eo and high-frequency €00 dielectric 
constants, and with the frequencies of the c axis IR- and Raman-active optical phonons j3| ■ 
The resulting Buckingham potential parameters, ionic shell charges {Y) and force constants 
{k) are presented in Table 7. To simulate different in-plane and out-of-plane bond lengths 
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TABLE VII: Potential parameters 
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147997.9 
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n nn 
U.UU 


Cul-Sr 


167925.4 


0.22873 


0.00 


Sr-Sr 


3680.5 


0.25880 


0.00 


Sr-Bi 


158000.0 


0.22873 


0.00 


Bi-Bi 


7100.0 


0.22873 


0.00 


Ca-Ca 


4500.0 


0.22873 


0.00 


0-0 


22764.0 


0.14900 


75.00 



Ion Y (I e I) k (eV -A-^) 

Ca 2.00 99.0 

Cul, Cu2 2.00 999999.0 

Bi 2.33334 27.2 

Sr 5.09995 146.7 

01 -3.25760 26.0 

04 -3.25760 76.0 

02, 03 -3.25760 100.0 
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TABLE VIII: A comparison between the experimental and calculated structures of Bi2212, Bi2223, 
and Bi2201 {I^/mmm (D^j^)). Experimental values of the fractional atomic coordinates and corre- 
sponding calculated values are given, with the experimental errors for Bi2212 structure shown in 
brackets. 





Bi2212 




Bi2223 




Bi2201 






z, exp. 


z, calc. 


z, exp. 


2;,calc. 


z, exp. 


z, calc. 


Ca 


0.00 


0.00 


0.0443 


0.0464 






Sr 


0.1097(5) 


0.1096 


0.1346 


0.1369 


0.0759 


0.0696 


Bi 


0.3022(3) 


0.3043 


0.292 


0.297 


0.3148 


0.3175 


Cul 


0.4456(9) 


0.4334 


0.4109 


0.4010 


0.00 


0.00 


Cu2 






0.00 


0.00 






01 


0.446(3) 


0.445 


0.411 


0.408 


0.00 


0.00 


02 


0.375(4) 


0.372 


0.353 


0.352 


0.405 


0.401 


03 


0.205(4) 


0.211 


0.213 


0.219 


0.194 


0.201 


04 






0.00 


0.00 







in the structures of Bi-based compounds, as well as anisotropy of oxygen polarizability, 
generally we have arrived at different pair potentials for the Cu-0, Sr-0, and Bi-0 in- 
plane and out-of-plane short-range interactions. It is particularly remarkable that the pair 
potentials for the Bi-02 and Bi-03 short-range interactions are significantly different, and 
this fact seems to be mainly due to very different Bi-0 bond lengths, i.e., in the structures 
of Bi-based compounds we have the Bi-02 bond length of 2.22 A, the in-plane Bi-03 bond 
length of 2.71 A, and the longest Bi-03 bond length, between weakly coupled Bi202 biplanes, 
of 2.97 A j^]. Thus elaborated, our shell model parameters allow us to describe all three 
Bi-based structures by the same set of short-range pair-potential parameters, listed in Table 



7. To compare the experimental and calculated structures we 



27 



present experimental values 



28l | and the corresponding 



of the fractional atomic coordinates obtained from Refs. 
calculated values in Table 8, where we also show the experimental errors for the Bi2212 
fractional coordinates in brackets. As far as our final fit is based on the Bi2212 crystal, 
the best agreement with the internal structural parameters is naturally obtained for this 
system. As can be seen from Table 8 these parameters reproduce also reasonably well the 
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Bi2223 and Bi2201 structures. Slightly varying the shell model parameters we can arrive at 
a better agreement with the particular structure chosen, although we would rather aim to 
obtain a self-consistent set of shell model parameters, which can account for the common 
generic features in all three Bi-based structures. Going ahead, we would like to mention here 
that this would not qualitatively affect the character of the eigenmodes resulting from our 
shell model calculations, although the eigenfrequencies and amplitudes of the eigenvectors 
of the c axis phonon modes can slightly change. The elastic properties calculated with these 
parameters do not exhibit unphysical (i.e., negative values), and the eigenvalues for the first 
three modes have approximately zero values in all three compounds which is indicative of a 
stable structure. We have also arrived at a good agreement with the experimental values for 
dielectric constants. All three compounds Bi2212, Bi2223, and Bi2201 have close values of c 
axis dielectric constants within the range eo — 9.5 -^ 12.0 and eoo — 4.9 -^ 5.8, as estimated 
from our experimental ei(a;) and €2(00) spectra. The corresponding calculated values fall 
in the range eo ~ 10.9 13.9 and — 4.8 5.6. Certainly, the obtained set of shell 
mode parameters is not unique. Nevertheless, given the fact that by using the same set 
of parameters we can reproduce the structures and main phonon bands as observed for all 
three compounds we obtain more confidence in the assignment for the phonon modes. 

B. Calculated T-point A2„ symmetry modes in Bi2201, Bi2212 and Bi2223 

In Figure 7 we show the eigenvector patterns corresponding to six eigenmodes of A2„ 
symmetry and one silent mode of B2U symmetry for the Bi2212, resulting from our lattice- 
dynamical calculations. Under each patch we list values of the eigenfrequencies and below 
we make an assignment based on a comparison with the experimental frequencies of the TO 
phonon bands observed in the c-polarized spectra of the Bi2212 single crystal (T = 150 K, 
see Figure 3 and Table 6). Mixing due to out-of- phase motion of ions of equal symmetry 
in inequivalent positions often does not allow the determination of a unique character for a 
given set of equal symmetry modes. The total number of the modes in each symmetry, of 
course, is unaffected by mixing. The character of the assigned modes is determined according 
to the dominant contribution from the atoms participating in the vibration and shown atop. 
The respective calculated components of the eigenvector amplitudes are given in Table 9. 
The lowest frequency mode calculated at 93 cm~^ (experimentally observed at 97 cm~^, see 
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Figure 3 and Ref. [5|) is related to motion of the heavy Bi atoms, which is strongly mixed 
with out-of-phase motion of Sr, Cul and Ca, with a small oxygen contribution from different 
crystallographic planes. The next mode of symmetry calculated at 195 cm~^ (168 cm~^) 
is predominantly due to motion of the Cul against Sr, with a small contribution from the 
Ca ions. The mode at 304 cm^^ (210 cm^^) is due to the vibrations of Ca and 01 
ions against Sr and 03. This mixed low-frequency mode contains noticeable contribution 
from the copper-plane 01 ions and bismuth-plane 03 ions to the c axis oscillations. The 
three high-energy vibrations are predominantly oxygen-related. The mode calculated at 343 
cm~^ (304 cm^^) is basically due to 03 and 01 vibrations against Sr and Ca atoms, with 
a predominant contribution from the bismuth-plane oxygen 03. The Cu-0 bond-bending 
character mode at 403 cm~^ (358 cm~^) involves significant contribution from the copper- 
plane 01 ions, vibrating against the Ca and 03 atoms. The highest frequency mode at 
643 cm~^ (583 cm~^) is due to 02 and 03 out-of-phase oscillations, and it has apparently 
apical character due to predominant contribution from the strontium-plane 02 oscillations 
along the c axis. Our assignment for the c-polarized phonon modes is confirmed by the 
oxygen isotope effect in the Bi2212 single crystal, investigated in this work (see Figure 
5). A noticeable red shift upon ^^O substitution is observed for the four c axis phonons 
at 210, 304, 358, and 583 cm~^; the estimated values of the shifts are 4.1, 4.4, 5.3, and 
5.5%, respectively. This result is consistent with our assignment that suggests that the 
three high-frequency modes have predominantly oxygen character. In addition, according to 
our assignment, the low-frequency mode at 210 cm~"^ also contains appreciable contribution 
from the copper-plane 01 and bismuth-plane 03 atoms, thus consistent with the noticeable 
red shift with the ^^O for ^^O substitution. The present theoretical shell model calculations 
thus reproduce the frequencies of the c axis eigenmodes quite well for the bilayer Bi2212 
single crystal. We would like to notice that although the eigenfrequencies are somewhat 
different the character of vibrations resulting from the present shell model calculations 
'or the Bi2212 system is quite consistent with the earlier calculations, made by Prade et al. 

In Figure 8 we show the eigenvector patterns corresponding to the eight eigenmodes 
of symmetry and two silent modes of B2U symmetry for the Bi2223, resulting from our 
lattice-dynamical calculations with the same set of the shell model parameters listed in Table 
7. The calculated components of the eigenvector amplitudes are summarized in Table 9. We 
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make an assignment of the TO phonon modes observed in the c axis elhpsometry spectra 
of the Bi2223 single crystal (T = 150 K, see Figure 3 and Table 6) based on a comparison 
with the present calculations. As can be seen from Figures 7 and 8, close similarity with 
the Bi2212 is observed for their eigenmode characters and eigenfrequency values due to 
the common generic features in the Bi2212 and Bi2223 structures. The only exception is 
that two new additional modes of symmetry are present for the trilayer Bi2223 system. 
Our calculations give one additional mode for Bi2223 at 145 cm~^, its character is defined 
by the main contribution from the Cu2 (from the inner Cu02 plane) and Sr out-of-phase 
vibrations, and assigned it with the new mode clearly discernible in the ei{uj) spectrum at 
127 cm~^ (see Figure 3(b)), which is also in agreement with the results of IR measurements 
done very recently on an oriented ceramics of Bi2223 We calculated four "light", 
predominantly oxygen-related, modes at the following frequencies 335, 378, 516, and 666 
cm~^. Like in the Bi2212, the highest-frequency mode at 666 cm~^ (581 cm^^) is associated 
with the out-of-phase 02 and 03 vibrations, the mode calculated at 335 cm~^ (304 cm~^) 
is predominantly due to the Bi-plane 03 oscillations, and the mode calculated at 378 cm~^ 
(359 cm~^) is the Cu-0 bond-bending character mode, corresponding to the out-of-phase 
motion of the Ca atoms and oxygens 01 in the outer CUO2 planes. We calculated one 
additional oxygen- related mode in the Bi2223 at 516 cm~^, and its character is defined by 
the out-of-phase vibrations of the Ca atoms and oxygens 04 in the inner CUO2 plane. We 
assign it with the new mode appearing in c axis the phonon spectrum of the Bi2223 at 402 
cm~^, as it is clearly observable in Figure 3, and also in agreement with the results of jiol . 

As we discussed in Section 4 some of the c-polarized IR phonon modes in the bilayer 
Bi2212 and trilayer Bi2223 compounds exhibit a strongly anomalous temperature depen- 
dence at Tc, accompanied with the development of a sizable absorption peak below Tc in 
the FIR range, as illustrated in Figure 6 . The two most pronounced phonon bands 

at 304 and 582 cm~^ appear at the same frequencies in the c axis ellipsometry spectra of the 
Bi2212 and Bi2223. According to our assignment they are associated with the predominant 
contribution from the BiO layer 03 oxygens and the apical 02 oxygens, respectively. It is 
evident from Figure 6 that the formation of the t-JPR in Bi2212 and Bi2223 is accompanied 
with the strongly anomalous temperature dependence of the phonon modes at 360, 402 (spe- 
cific for trilayer Bi2223) and 582 cm~^ denoted by A, B, C. An explanation of the anomalous 
changes of the phonon modes at 360 and 402 cm~^ upon entering the superconducting state 
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has been recently proposed based on the so-called Josephson superlattice model, assuming 
that these bands belong to the vibration of the planar oxygens jsl [t]. This assumption is 
supported by our assignment which suggests that these modes are bond-bending character 
modes which involve a significant contribution from the outer copper-plane 01 and inner 
copper-plane 04 ions, respectively. 

In Figure 9 we present the eigenvector pattern corresponding to the five eigenmodes 
of symmetry for Bi2201, the calculated components of the eigenvector amplitudes are 
listed in Table 9. We have calculated two symmetry modes at low frequencies. The 
calculated lowest frequency mode is related to the motion of the heavy Bi atoms, which is 
strongly mixed with out-of-phase motion of Sr and Cul. The next mode of symmetry is 
predominantly due to the motion of Cul against Sr. The three phonon features are clearly 
observed in the c axis ellipsometry spectra of the Bi2201 at low frequencies at around 109, 
165, and 209 cm~^ (see Figure 4 and Table 6); in Section 4 we have already discussed 
the possible reasons for the appearance of one "extra" phonon feature. We have made an 
assignment of the modes at low frequencies by analogy with the vibration characters at 
the corresponding frequencies in the B2212 and Bi2223 crystals. The three high-energy 
vibrations are predominantly oxygen-related. The highest frequency mode at 596 cm~^ 
(586 cm^^) has apparently apical character due to the predominant contribution from the 
strontium-plane 02 oscillations along the c axis. The modes calculated at 293 cm~^ and 353 
cm~^ are basically due to bismuth-plane 03 and copper-plane 01 mixed in-phase and out- 
of-phase vibrations, respectively, and they correlate well with the experimental frequencies 
at 301 and 385 cm~^. In the Bi2201 the CUO2 planar oxygen band located at 385 cm~^ has 
much larger intensity than in the Bi2212. This is very probable due to the fact that in the 
Bi2212 the CUO2 planar oxygen mode is, already in the normal state, strongly renormalized 
by its coupling to interlayer electronic excitations. 



C. Calculated T-point Aig symmetry modes in Bi2201, Bi2212 and Bi2223 

In this Section we just briefiy discuss the results of our shell model calculations for the 
Raman-active modes in the three Bi-based compounds and suggest the assignment for the 
phonon modes observed in the experimental spectra. It has been generally found that a larger 
number of phonon modes occur in the Raman spectra of the Bi-based compounds than would 
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TABLE IX: Calculated components of eigenvectors of A-2u modes in Bi2201 (Bi2Sr2Cu06), Bi2212 
(Bi2Sr2CaCu208) and Bi2223 (Bi2Sr2Ca2Cu30io) crystals in tetragonal I4/mmm (D^^) space 
group. 



Bi2212 






Cul(z) 




Ca(z) 




02{z) 


03(2;) 




93 CTTl^^ 


-0.46 


0.27 


0.33 




0.29 


0.16 
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-0.11 




195 cm~^ 
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0.00 




304 cm-^ 
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0.18 




343 cm~^ 
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-0.18 
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Bi2223 
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Cu2(z) 


Ca(z) 


Oliz) 
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04(z) 


88 cm~^ 


-0.45 


0.11 


0.14 


0.57 


0.21 


0.08 


-0.07 


-0.11 


0.17 


1 zL,^ fin 


-0.23 


0.33 


0.22 


-0.68 


0.09 


0.10 


0.00 


-0.01 


-0.09 


211 cmT^ 


0.01 


0.27 


-0.59 


-0.06 


0.23 


0.06 


-0.03 


0.00 


0.09 


252 cm"^ 


-0.05 


0.42 


-0.05 


0.29 


-0.49 


-0.07 


0.02 


0.03 


-0.16 


335 cm~^ 


-0.09 


-0.15 


-0.07 


-0.01 


-0.21 


0.29 


0.10 


0.46 


0.16 


378 cm~^ 


0.10 


-0.09 


-0.05 


0.09 


-0.13 


0.34 


-0.14 


-0.41 


0.04 


516 cwT^ 


0.05 


0.05 


0.03 


-0.21 


-0.23 


-0.05 


-0.08 


-0.22 


0.60 


666 cm~^ 


0.08 


0.03 


0.05 


0.00 


0.01 


-0.01 


-0.66 


0.22 


0.02 


Bi2201 


Bi(z) 


Sr(z) 


Cul(z) 






Ol{z) 


02(z) 


03(z) 




77 cm"-^ 


-0.31 


0.24 


0.77 






0.10 


-0.11 


-0.17 




175 cmT^ 


-0.21 


0.51 


-0.56 






0.17 


0.02 


-0.06 




293 cm-^ 


-0.24 


-0.09 


0.08 






0.30 


0.15 


0.56 




353 cm~^ 


0.10 


-0.23 


-0.08 






0.59 


-0.10 


-0.27 




596 cm~^ 


0.09 


0.06 


-0.05 






0.00 


-0.66 


0.23 





be expected from the simple tetragonal M/mmm symmetry consideration (given in Tables 
1 and 2). As a result different groups have reported rather different and then contradictory 
assignments. We have already referred to these reasons in Section 2. There appears to be 
a consensus on the origin of the strong Bi^ phonon in Raman spectra of Bi2212. From 
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the Big-symmetry selection rules the peak at about 285 cm~^, which occurs only in the 
z{xy)z spectrum can be assigned to the 0(1) Bi^ phonon characterized by out-of-phase 
motion of the oxygen atoms in the CUO2 plane. This is not the case, however, for two of 
the strongest Ai^ modes that appear in the Raman spectrum at higher frequencies around 
460 and 620 cm~^. The origin of the high-frequency vibrations in the Raman spectrum of 
Bi2212 remains controversial despite many investigations by different groups. Most authors 
have assigned the 460 cm~^ mode to the vibrations of the apical oxygen atoms located in the 
SrO planes and the 620 cm~^ phonon to vibrations of the oxygen atoms located in the BiO 



olanes BMtal. However, reeeutlv ,t was argued that the assignment should be reversed 
|48|,|49||. It was shown that these two bands exhibit different polarization dependencies for 
the incident and scattered vectors, indicating that the degree of mode mixing is small, and 
the 620 cm~^ band is strongly c -polarized like apex oxygen vibrations characteristic for 
most high-Tc superconductors [48]. Moreover, the 460 cm~^ band has been found to soften 
at a significantly faster rate than the 620 cm~^ band with ^^O oxygen isotope substitution, 
in agreement with the more liable oxygen atoms from the BiO layer [4^. Therefore, in 
fitting our shell model parameters for the Bi2212 crystal we have implied an assignment for 
the mode observed at around 285 cm^^ as the Bi^ mode, whereas for the modes observed 
at around 460 and 620 cm"^ as the Aig modes associated with the 03 (from the BiO plane) 
and 02 (from the SrO plane) predominant contributions, respectively. 

In Figure 10 we show our suggested eigenvector pattern corresponding to the six eigen- 
modes of Ai^ symmetry and one Big symmetry for the Bi2212 crystal. We list values of the 
eigenfrequencies and below we make an assignment in agreement with the experimental Ra- 
man study of Bi2Sr2Cai_^Yj.Cu208+5 (x = 1) single crystals [4Sj. A number of phonons 
observable below 400 cm~^ in the Raman spectra in the Aig (or Ag) scattering configura- 
tions are considered by Kakihana et al. as " disorder- induced" modes 3|. In addition, as a 
weak satellite feature at ~ 657 cm~^ near the strong phonon band at 620 cm~^ disappears 
in modulation-free Bi2_xPbxSr2Cai_^Cu208+5 crystals it is possibly induced by the incom- 
mensurate superstructural modulation j^. Another satellite feature at ~ 458 cm~^ near 
the pronounced phonon band at 460 cm~^ is possibly associated with extra oxygen atoms 
in the BiO layers, and also referred to " disorder- induced" scattering. 

In Figure 11 we show the eigenvector pattern corresponding to the seven eigenmodes 
of Aig symmetry and one Big symmetry for the Bi2223 crystal, calculated using the same 
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basic set of the shell model parameters. From Figures 10 and 11 we can see that due to 
the common generic features in the Bi2212 and Bi2223 structures the vibration characters 
are similar and the eigenfrequency values are close for the corresponding eigenmodes. This 
is consistent with the recently published results on Raman measurements showing that 
frequencies of the major phonon bands are almost identical in Bi2212 and Bi2223 single 
crystals [s^]. The exception is an additional mode of the Aig symmetry present in the 
trilayer Bi2223 crystal. Comparing the bilayer and trilayer systems (Figures 10 and 11) we 
can conclude that instead of one Ai^ mode associated with the 01 vibrations in the Cu02 
plane in the Bi2212 structure, calculated at 370 cm~^, in the Bi2223 structure we have 
two Aig modes associated with the Ca and 01 (from the outer Cu02 planes) in-phase and 
out-of-phase vibrations, calculated at 316 and 391 cm~^, respectively. 

In Figure 12 we show the eigenvector pattern corresponding to the four eigenmodes of Ai^ 
symmetry for the M/mmm Bi2201, resulting from our shell model calculations. As usual 
many more phonon modes are observed in the experimental Raman spectra implying that the 
simple tetragonal approximation is no more applicable for the interpretation of the Raman 
phonon spectra of the Bi2201. In Tables 1 and 3 we show that due to the orthorhombic 
distortions the 01 atoms in the Cu02 plane become "Raman-active", and, as a result, a 
number of new Raman-active modes of A^, Bi^, 2B2g, and 2B3g character will appear in the 
Amaa Bi2201. We suggest that a more appropriate approximation for the interpretation of 
the F-point Raman-active phonons in this compound will be the orthorhombic Amaa space 
group. However, the assignment for the two high-frequency oxygen modes (as we show in 
Figure 12) is the same as for the Bi2212 and Bi2223 crystals, and the calculated frequencies 
agree quite well for these three compounds. 

Resonant Raman scattering below Tc has been discovered in Bi- and Hg-based high Tc 
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53|. Strong superconductivity- induced 



compounds with three or four Cu02 layers 
enhancement of the phonons at 260 and 390 cm~^, together with the development of a 
broad band peaking at around 580 cm~^, has been observed in the recent resonant Raman 
scattering study on the Bi2223 single crystal with Tc = 109 K j^^. Our measurements of 
resonant Raman scattering (E^xc = 2.18 eV) done on our Bi2223 single crystal with lower 
Tc = 107 K reproduce the results of Limonov et al. In Figure 13(a) we show the Raman 
spectra at two different temperatures, i.e. somewhat above and well below Tc. Figure 
13(b) shows the difference between these Raman spectra. Our data clearly demonstrate 
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the superconductivity-induced enhancement of the phonon bands at 260 and 390 cm~^ and 
a broad band peaking around 550 cm~^. There are two ahernative explanations of these 
superconductivity-induced anomahes in multilayer high- compounds; one considers the 
broad excitation band as a "pair-breaking" peak at 2 A [52], and the another refer it to a 
Raman-active c axis plasmon IssJ. A more detailed study of the doping effect on the res- 
onance peak position is required to clarify its origin. Nevertheless, dicussing the phonon 
anomalies observed in trilayer Bi2223 we would like to notice that both approaches consider 
the superconductivity-induced anomalies as due to strong coupling of the intratrilayer elec- 
tronic excitation with phonons which involve calcium and planar oxygen ions. Experimental 
resonant Raman scattering studies bring evidence in favour of our assignment suggested 
for the trilayer Bi2223 crystal (Figure 11). The superconductivity-enhanced phonon bands 
at 260 and 390 cm~^ in the Bi2223 can be ascribed to the Ai^ in-phase and out-of-phase 
vibrations of the Ca atoms with oxygens in the outer Cu02 planes, in agreement with our 
assignment for the Ai^, modes calculated at 316 and 391 cm~^. 



VI. CONCLUSIONS 



In summary, we present experimental data on single crystals as obtained by spectral el- 
lipsometry and make the assignment of the IR c axis polarized phonon modes in single, bi-, 
and trilayer bismuth compounds (Bi2201, Bi2212, and Bi2223) by comparing them with the 
results of the shell model calculations. The c axis IR phonon spectra look very similar for 
the bilayer and trilayer Bi-based crystals, supporting experimental evidence for the generic 
structural features of these compounds. We apply a tetragonal lA/mmm approximation in 
the interpretation of the c axis IR spectra of the Bi-based compounds, which are known to 
have pseudotetragonal structure, and conclude that this approximation provides a reasonable 
description of the main features in the phonon spectra. In the framework of this approx- 
imation, we have identified " disorder- induced" modes due to the orthorhombic distortions 
associated with the incommensurate superstructural modulations and extra-stoichiometric 
oxygen. Comparing with the bilayer Bi-system, we identify two additional c axis phonon 
modes in the ellipsometric spectra of the trilayer system, one at 127 cm~^, among the "heavy" 
modes, and the another at 402 cm"^, among the "light" modes, which according to our as- 
signment involve significant contribution of the copper and oxygen atoms in the inner Cu02 
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plane, respectively. The recent experimental observation of a Josephson-plasma-resonance 
and phonon anomalies in Bi2223 7[ provide, indeed, a strong support for the Josephson- 
superlattice model, which consistently explains, the phonon anomalies and the formation 
of a broad absorption band below Tc using the components of the eigenvector amplitudes 
for the relevant oxygen modes at 358 and 402 cm^^. The two modes involving 03 (BiO 
plane) and 01 (Cu02 plane) atoms exhibit substantial frequency changes in the single layer 
system, and appear at 294 and 389 cm~^, respectively. For the single layer system the mode 
at 389 cm^^ increases dramatically in intensity. This evident difference between single and 
multilayer compounds in the behavior of this particular phonon mode reflects the fact that 
in multilayer compounds the oxygen bond-bending phonon mode is strongly-coupled with 
the electronic excitation due to the charge density oscillations between the closely spaced 
Cu02 planes. In addition, we also show the eigenvector patterns and calculated frequencies 
for the Raman-active phonon modes of Ai^ symmetry in these three compounds according to 
the results of the present shell model calculations. The supeconductivity-induced resonant 
Raman scattering for the two phonon modes at 260 and 390 cm~^ very recently observed in 
the trilayer Bi2223 single crystal [3| is well consistent with the suggested assignment. 
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Bi^Sr^Ca^CuP; 
14/mmm 




FIG. 1: Crystallographic cells of Bi2Sr2Can_iCun02n+4 with n = 1, 2 and 3. 
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FIG. 2: Notations for copper and oxygen atoms in primitive unit cells of single, bi-, and trilayer 
Bi-based cuprates. 
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FIG. 3: Real parts of the c axis (a) optical conductivity ai{u) and (b) dielectric function ei{u) of 
Bi2212 and Bi2223 single crystals with Tc = 91 K and Tc = 107 K at T = 150 K. 
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FIG. 4: C axis (a) optical conductivity 0-1(0;) and (b) dielectric function ei{uj) of Bi2201 single 
crystal at T = 150 K. 
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FIG. 5: Oxygen isotope effeet in the c-polarized response of (a) (7i(u)) and (b) €i(a;) of Bi2212 with 
Tc = 91 K, measured at 150 K. 
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FIG. 6: Real part aUio) of the FIR c axis conductivity (a) of Bi2212 with = 91 K for T = 
100 K and T = 10 K 6] and (b) of Bi2223 with = 102 K for T = 120 K and T = 10 K |2|. 
(c) Superconductivity-induced changes of the FIR conductivity, Aai = (Ti{T « TcUj) — (Ti{T > 
Tc, io), for Bi2212 and Bi2223 derived from the §§ta shown in the top panels. The most pronounced 
phonon anomalies are denoted by A, B, and C. 



Infrared Aj, (Bi221 2 14/mnnm) 

Bi':CulCaSr Sr:Cul' Co : Sr' 0301 01':Ca03 02:03' Silent B, 




93 195 304 343 403 643 230 



(97) (168) (210) (304) (358) (583) 

FIG. 7: Normal frequencies (in cm^^) of six symmetry modes and one "silent" mode in 
Bi2212 [lA/mmm) and schematic representation of their eigenvectors; respective experimental c 
axis TO phonon frequencies are given in brackets (see Figure 3 and Table 6). 
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Infrared A^, (Bi2223 14/mmm] 
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FIG. 8: Normal frequencies and eigenvectors of eight symmetry modes and two "silent" 
modes in Bi2223 (M/mmm); respective experimental c axis TO phonon frequencies are given in 
brackets (see Figure 3 and Table 6). 
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Infrared A^, (Bi2201 14/mmm] 
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FIG. 9: Normal frequencies and eigenvectors of five symmetry modes in Bi2201 (M/mmm); 
experimental c axis TO phonon frequencies are given in brackets (see Figure 4 and Table 6 ). 



Ramon A,g (Bi221 2 14/mmnn) 

BiSr CulSr SrCul' oi:Sr' 03 02 B, 




112 149 208 370 461 630 273 

(59) (117) (145) (409) (463) (627) (287) 



FIG. 10: Normal frequencies and eigenvectors of six Ai^ symmetr y rn odes and one Bi^ mode in 
Bi2212 {lA/mmm); experimental frequencies are given in brackets |48l [. 
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Raman A,g (Bi2223 14/mmm) 

SrCulBi Bi:Sr' Cul:Sr' Ca01:Sr' Ca:01' 03 02 B, 




103 125 214 316 391 446 646 237 

(260) (390) (463) (627) (287) 



FIG. 11: Normal frequencies and eigenvectors of seven Aig symmetry modes and one Bi^ mode in 
.3i2223 (lA/mmm); some of experimental frequencies are given in brackets (see Figure 13 and Ref. 

3)- 



Raman A,g (Bi2201 14/mmm) 
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FIG^ 12: No...al .e,ue„cies and ei„s of ta. A, s™.e.. .ode. i. Bi2201 (.4/™,; 

some of the experimental frequencies are given in brackets |48|. 
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FIG. 13: Raman scattering {Eexc = 2.18 eV) in Bi2223 (a) measured slightly above Tc = 107 K 
and well below Tc and (b) resonant enhancement of the phonon bands at 260 and 390 cm~^ and 
superconductivity-induced broad band at 550 cm^^. 
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